Self-intersection local time of planar Brownian 
motion based on a strong approximation by 
random walks 



Tamas Szabados* 
Budapest University of Technology and Economics 

March 3, 2011 



Abstract 

The main purpose of this work is to define planar self-intersection local 
time by an alternative approach which is based on an almost sure pathwise 
approximation of planar Brownian motion by simple, symmetric random 
walks. As a result, Brownian self-intersection local time is obtained as an 
almost sure limit of local averages of simple random walk self-intersection 
local times. An important tool is a discrete version of the Tanaka-Rosen- 
Yor formula; the continuous version of the formula is obtained as an almost 
sure limit of the discrete version. The author hopes that this approach to 
self-intersection local time is more transparent and elementary than other 
existing ones. 



1 Introduction 

Let (W(t))t>o be planar Brownian motion (BM). Formally, its self-intersection 
local time at the point x G K 2 up to time t is 

a(t,x)= ( [ S(W{v) -W{u) - x)dudv, 
Jo Jo 

where S is the Dirac measure at zero. There exist several methods in the litera- 
ture to make this definition rigorous. One natural approach, which is the topic 
of the present work, is to define a(t, x) as an almost sure limit (when m — > oo) 
of local averages of self- intersection local times a m (t, x) of a nested sequence of 
simple, symmetric planar random walks (B m (t)). (See the next section for the 
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definition of B m .) For these imbedded random walks (RW's), self- intersection 
local time can be defined by elementary counting: 

a m (t,x) := 2" 2m #{(i,i) : < i < j < t2 2m , B m { 3 T 2m ) - B m {i2~ 2m ) = x}. 

Then, as it will be seen in Theorem [21 the following almost sure limit gives the 
Brownian self-intersection local time: 

a(t,y)= lim lim -!- V a m (i,x)2" 2m (y ? 0), 

(5->0+ nwoo 7T0 * — ' 

xe2-™1 2 C\B s (y) 

where Bg(y) denotes the disc centered at y with radius 5. 

The author hopes that this approach to self-intersection local time is more 
elementary and more advantageous from a pedagogical point of view than other 
existing ones. This method is a special case of a strong invariance principle for 
self-intersection local time. It should be mentioned that using different methods, 
strong invariance was shown earlier by Cadre [5] , and for general random walks 
by Bass and Rosen pQ . The method applied in this paper is based on a Tanaka- 
like formula, first introduced by Rosen [13] and then generalized by Yor [2D] . 
More exactly, a discrete version of the Tanaka-Rosen-Yor formula for random 
walks is given below whose almost sure limit is the continuous version of the 
formula. 

2 Preliminaries 

A basic tool of the present paper is an elementary construction of Brownian 
motion. The specific construction used in the sequel, taken from |T7], is based on 
a nested sequence of simple, symmetric random walks that uniformly converges 
to the Wiener process (=BM) on bounded intervals with probability 1. This 
will be called "twist and shrink" construction. This method is a modification 
of the one given by Frank Knight in 1962 [7]. 

We summarize the major steps of the "twist and shrink" construction here. 
We start with an infinite matrix of independent and identically distributed ran- 
dom variables X m (k), P {X m (k) = ±1} = | (m > 0, k > 1), defined on the same 
complete probability space (SI, J-, P). (All stochastic processes in the sequel will 
be defined on this probability space.) Each row of this matrix is a basis of an ap- 
proximation of the Wiener process with a dyadic step size At — 2~ 2m in time and 
a corresponding step size Ax = 2~' m in space. Thus we start with a sequence 
of independent simple, symmetric RW's S m (0) = 0, S m (n) = Y^k=i^m[k) 
(«>!)• 

The second step of the construction is twisting. From the independent RW's 
we want to create dependent ones so that after shrinking temporal and spatial 
step sizes, each consecutive RW becomes a refinement of the previous one. Since 
the spatial unit will be halved at each consecutive row, we define stopping times 
by T m (0) = 0, and for k > 0, 

T m (k + 1) = min{n : n > T m (k), \S m {n) - S m {T m (k))\ =2} (to > 1) 

These are the random time instants when a RW visits even integers, different 
from the previous one. After shrinking the spatial unit by half, a suitable modi- 
fication of this RW will visit the same integers in the same order as the previous 
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RW. We operate here on each point u> e SI of the sample space separately, i.e. 
we fix a sample path of each RW. We define twisted RWs S m recursively for 
k = 1,2, ... using S m -i, starting with So(n) = So(n) (n > 0) and S m (0) = 
for any m > 0. With each fixed m we proceed for k = 0,1,2,... successively, 
and for every n in the corresponding bridge, T m (k) < n < T m (k + 1). Any 
bridge is flipped if its sign differs from the desired: 



X m (n) = 



X m {n) if S m (T m (k + 1)) - S m {T m (k)) = 2A > m _ 1 (fc + 1), 
-X m (n) otherwise, 



and then S m (n) = S m (n — 1) + X m (n). Then S m (n) (n > 0) is still a simple 
symmetric RW |T7l Lemma 1]. The twisted RWs have the desired refinement 
property: 

S m +i{T m +i{k)) = 2S m {k) (m > 0, k > 0). 

The third step of the RW construction is shrinking. The sample paths of 
S m (n) (n > 0) can be extended to continuous functions by linear interpolation, 
this way one gets S m (t) (t > 0) for real t. The mth "twist and shrink" RW is 
defined by 

B m (t)=2- m S m (t2 2m ). 
Then the refinement property takes the form 

B m+ i (T m+1 (fc)2- 2(m+1) ) = B m (k2- 2m ) (m > 0,k > 0). (1) 

Note that a refinement takes the same dyadic values in the same order as the 
previous shrunken walk, but there is a time lag in general: 

T m+ i(fc)2~ 2(m+1) - k2- 2m ^ 0. (2) 

It is clear that this construction is especially useful for local times, since a re- 
finement approximates the local time of the previous walk, with a geometrically 
distributed random number of visits with half-length steps, cf. |18j . 

Now let me recall some important facts from [T7] and [TH] about the "twist 
and shrink" construction that will be used in the sequel. 

Theorem A. On bounded intervals the sequence (B m ) almost surely uniformly 
converges as m — > oo and the limit process is Brownian motion W . For any 
C > 1, and for any K > and m > 1 such that K2 2m > N{C), we have 

sup \W{t) - B m {t)\ > 27C#jF(log„ K)imh- 

0<t<K 

where := K V 1 and log^ K := (log if) V 1. 

(N(C) here and in the sequel denotes a large enough integer depending on 
C, whose value can be different at each occasion.) 

Conversely, with a given Wiener process W, one can define the stopping 
times which yield the Skorohod embedded RWs B m (k2~ 2m ) into W. For every 
m > let s m (0) = and 

s m (k + l) = inf{s: s> s m (k),\W(s) - W{s m {k))\ = 2~ m } (k > 0). (3) 
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With these stopping times the embedded dyadic walks by definition are 

B m {k2- 2m ) = W(s m (k)) (m > 0, k > 0). (4) 

This definition of B m can be extended to any real t > by pathwise linear 
interpolation. 

If a Wiener process is built by the "twist and shrink" construction described 
above using a sequence (B m ) of nested RWs and then one constructs the Sko- 
rohod embedded RWs (B m ), it is natural to ask about their relationship. The 
next theorem shows that they are asymptotically equivalent. In general, roughly 
saying, (B m ) is more useful when someone wants to generate stochastic processes 
from scratch, while (B m ) is more advantageous when someone needs discrete 
approximations of given processes. 

Theorem B. For any C > 1, and for any K > and m > 1 such that 
K2 2m > N{C) we have 

sup \W(t)~B rn (t)\ > 27CK}{\og,K)imi2-^ 

0<t<K 

< Y^c(^ 2m r c - 

Apply the "twist and shrink" construction <i-times independently, to obtain a 
d-dimensional Brownian motion W = (W 1 , . . . , W d ) (vector components will be 
denoted by superscripts), the corresponding Skorohod-embedded RWs B m = 
(Bl n , . . . , B^), and stopping times (s I 1 n (n), . . . , s%): 

Bl(k2- 2m ) = W j {sl{k)) (m > 0, k > 0,j = 1, . . . ,d). 

Please note that in this paper d-dimensional random walks are defined as a 
vector of d independent one-dimensional random walks. This means that the 
coordinate axes of a usual (i-dimensional random walk are rotated and the length 
of a step is multiplied by \fd. 

Then Theorem IB1 and Borel-Canteli lemma imply 

Corollary 1. Taking a d- dimensional Brownian motion W and Skorohod em- 
bedded RWs B m , for any K > and m > 1 one has 



sup 

0<t<K 



\W{t)-B m {t)\ = O ((logn)i n"i) = o(m*2-^\ a.s., (5) 



where n — K2 2m denotes the number of vertices of an imbedded random walk 
B m over the time interval [0, K] . 

As it was mentioned above, the above approach to Brownian motion is es- 
pecially suitable to give an elementary definition of Brownian local time as an 
a.s. limit of RW local times, cf. Q~8]. In fact, this was the main motivation to 
find a similar definition of planar self-intersection local time as well. 

The idea that random walk approximations can be applied to obtain results 
about Brownian local time goes back to Knight [8], who proved the celebrated 
Ray-Knight theory this way. Revesz [TT] and Csaki & Revesz [3] were the first 
to prove strong invariance for local times. 
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Let £ m (0, x) := and 

£ m (k, x) := #{j : < j < k, S m (j) = x}. (6) 
Define local times of the "twist and shrink" RW B m or imbedded RW's B m by 
£ m (t,x) :=2- m £ m (t2 2m ,x2 m ). (7) 

Then [TS] shows 

Theorem C. On any sfrzp [0,-ftT] x M, 

lim £ m (t, x) = a;) a.s., 

m— >oo 

uniformly in (t,x), where C(t,x) is the local time of BM. Hence this automati- 
cally gives a version of Brownian local time which is continuous in (t, x) . 
Moreover, for all K > and m > 1 one has 



sup 

{t,x)e{0,K]xl 



£(t,x) - C m (t,x) =C>((logn)4n M 



where n = K2 Z 



Similar statements hold for £ m (t,x) computed from Skorohod embedded 
RW's B m as well. Interestingly, the rate of convergence is the same for local 
time as for the approximation of Brownian motion with the "twist and shrink" 
construction. While this rate is much weaker than the optimal Komlos-Major- 
Tusnady rate (logn) n~2 in the case of BM, it just slightly differs from the 
optimal (logn)a (log log n) i n~i in the case of local time. 



3 A discrete Ito's formula 

It is interesting that one can give discrete versions of Ito's formula and of Ito 
Tanaka-Meyer formula, which are purely algebraic identities, not assigning any 
probabilities to the terms. Despite this, the usual Ito's formula follows fairly 
easily in a proper probability setting. 

Discrete Ito formulas are not new. Apparently, the first such formula was 
given by Kudzma in 1982 [5]. The elementary algebraic approach used in the 
present paper is different from that; it was introduced by the author in 1989 

ESI- 

First we need definitions of discrete line integrals and conservative vector 
fields on a grid. Fix an initial point a £ M. d and step-size (mesh) h > 0. 
Consider the grid G(a, h) := a + hZ d , and let / = (J 1 , . . . , f d ) : Q(a, h) -> R d 
be a vector field on this grid. (Coordinates of a vector will always be denoted 
by superscripts.) Take an arbitrary broken line (a discrete path) 7 that goes 
through finitely many (not necessarily distinct) oriented edges between adjoining 
vertices of the grid. A typical such edge is [x, x + fxhej], where x € Q(a,h), ej 
(1 < j < d) is a coordinate unit vector and fi = ±1. (The order of the two 
vertices is important!) A discrete path 7 is a formal sum of such oriented edges 
(that is, a 1-chain): 

n 

7 = ^""^ [x r . x r + [i r hej r \ (1 < j r < d). 

r=l 
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By definition, the corresponding discrete path integral (or trapezoidal sum) 
of / over 7 is defined as 



h 

T 1 (/)7o) h: = 2^2^ (f^^r) + f 3r {Xr+ Vrhe jr )) . 
r=l 

Here the symbol 7q refers to the unit tangents of 7 along edges, and (/, 7o) 
denotes dot product. When it will be convenient, T 7 (f(x), "f'o(x)) h will be 
written to show the dummy variable of the summation. 

The ordinary path (line) integral of a vector field / over a path 7 will be 
denoted by J (/, 7q) ds or by f f(x) ■ da;, where ds refers to the length element 
of 7. 

The above definition of discrete path integral shows that the orientation of 
an edge is defined by the order of its two vertices: it is positive if the edge goes 
increasingly in a coordinate and negative in the opposite case. If 7 = 0, we 
define T 7 (f,y o )h= 0. 

A vector field / is called discrete conservative on the grid G(a, h) if for any 
b, c G Q{a,h) and for any discrete path 7 going from b to c through edges of 
adjoining neighbor vertices of Q(a, h), the discrete path integral does not depend 
on the path 7, it depends only on the initial point b and endpoint c. In this 
case the notation T h c (/, Jq) h will be used for the trapezoidal sum. Clearly, / is 
discrete conservative if and only if T§ (/, j' Q ) h — whenever b = c, that is, the 
path 7 is closed. 

When / is discrete conservative, one can define a discrete potential g : 
G(a, h) — > R by the formula g(x) :— (/, 7g) h. Then for any points b and c 
in the grid, and for any discrete path 7 connecting them, one has T 7 (/, 7q) h = 
9(c) ~9(b). 

The following discrete Itd's formula (which is a simple algebraic identity) al- 
ready appeared in [TH Section 5] in the two-dimensional case. It is based on the 
principle that though our random walk is "diagonal" , constructed from indepen- 
dent one-dimensional random walks, the discrete integrals below go "parallel to 
the coordinate axes" , as in a standard continuous Ito's formula. Also, the main 
object in our discrete formula is the "integrand" / in the stochastic sum, which 
corresponds to the gradient of a scalar field in a standard continuous Ito's for- 
mula. That may explain why we suppose that / be discrete conservative. These 
methods make it convenient to deduce important continuous formulae from the 
discrete ones. 



Lemma 1. Take a € Mr, step h > 0, and a discrete conservative time- dependent 
vector field f = (/ ,•••>/ ) : h 2 Z + x Q(a,h) — > R d . Consider a sequence 
X r = (X^,...,Xf) (r > 1), where X-j, = ±1. Define partial sums Sq = a, 
S n = a + h(Xi + ■ • ■ + X n ) (n > 1 ) and discrete time instants t r — rh 2 (0 < 
r < n). Assume that the steps of (S n ) are performed in time steps h 2 . Then 
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the following equalities hold: 

Tt So (f(tn,x), 1 ' (x))h 
n 

= E T "= So ({/(*r, X) - /(t r -l, x)} , i {x)) h 



r=l 



" * P (tr-i, Sr-i + Eti 1 ^ei) + /' (t r -i, 5 r _i + E-=i ^e,) 

+E E 5 



(8) 



r=l j = l 

(discrete Stratonovich formula) . Alternatively, 
T ts {fit n ,x),i Q {x))h 

n 

= Y. T t So ({f(tr,x) - f(t r - U x)} ,7o(*)) >» 
r=l 

n d / j-1 \ 

+ E E ^ *r-l. Sr-1 + E ^ 
r=l j=l V i=l / 

1 " * P (tr-1, Sr-1 + Ei=l - /' (t r -l, S r -1 + E£i hX* . , 

+ 2 E E — 

r— 1 ^ = 1 ' 

(9) 

(discrete Ito's formula). 
Proof. Algebraically, 

T x=S (f(tr,x),7o(x)) h - Tfl- (/(tr-1,*), Y (*)) ft 

= T^ So (f(t r ,x), 7 ' (x)) h - T x t So (/(ir-i.aO.YoCr)) ft 

+ r x t So (/(t r -l, X), 7 ' (x))h-T^- S l (f{tr-l,x),-/ Q (x))h. 

Using the assumption that / is discrete conservative, we get that 

Tt So (/(tr-l,!),^!))/!"^; (f(t r - 1 ,x),T , Q (x))h 

= T^ Sri (f(t r . u x),y (x))h 

and 

T^ Sr .Af(tr-i^)n'o(x))h (10) 



* P (*r-l, Sr-1 + ELl ^i) + / j (tr-1, Sr-1 + ELl 
E 9 hX r 



3 = 1 



= E ^ f*r-L S r-1 + E ^ e *) ^ 

j=i V »=i / 



1 * / j (tr-1, Sr-1 + ELl - P (tr-1, Sr-1 + ELl ^i) 

+ -> ^ : ^ '-h 2 . 

2 hXl 
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The first equality follows from the fact that if X£ — 1, one has a positively 
oriented edge, while if X-j, — — 1, one has a negatively oriented edge in the 
trapezoidal sum. Then the second equality follows since 1 jX 3 T = X d r . Summing 
up for r = 1, . . . , n, the sum on the left telescopes, and from the two equalities 
one obtains the two formulae, respectively. □ 



One can introduce partial local times Ci(t n ,x) (n > 0) of the series (S n ) 
with spatial step h > 0, time step h 2 , and fi G {1, —l} d : C^(0, x) := and 

C%(t n ,x) := h #{j : < j < n, Sj ; = x, S J+1 = x + hfi}, (11) 

where n > 1 and x £ Q(a,h). The (total) local time is 

C h (t n ,x):= tfl(t n ,x) = h#{j:0<j<n,S j =x}. (12) 

ne{i,-i} d 

Here our convention differs from the usual one: time is counted, but time n 
is not. The reason is that this better fits the discrete formula below. 

Lemma 2. With the same assumptions as above in Lemma\]\ except that the 
vector field f does not depend on time, f : Q(a,h) — > M. d , one also has 



T^ So (f(x),y (x))h (13) 



->s n 

n d / j—1 



r=l j=l \ i=l J 

+ \ E E 

KGa+hZ d fie{l,-l} d 

x E M j {/ j (x + E ^) - / J ( x + E 

(discrete Ito-Tanaka-Meyer formula). 

Proof. Continuing (|10[) in the proof of the previous lemma, 
= E f^-i + E WC? 

j=l V i=l / 

+ 2 E E ftl {£v_i=a!,£v=a:+h/i} 

E /'' { /' (•'• • E ''/'''-) - f j [■'■ + E ''i 1 '' • 

3=1 I V i=l / V i=l 



X 



Again, summing up for r — 1, . . . , n, the sum on the left telescopes, and on the 
right one obtains the asserted formula. □ 
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4 Constructing discrete conservative vector 
fields on a planar grid 

For the sake of simplicity, from now on only the planar case (d = 2) will be 
discussed, as this is the case that will be used in the sequel. The problem 
that we consider in this section is that given a differentiable scalar field g in 
the plane, its gradient Vg is not a discrete conservative vector field on a grid 
G(a, h) = a + hZ 2 in general. We want to modify Vg so that the resulting vector 
field f be discrete conservative on the grid, but still do not differ much from 

Let us call any E := [x,x + he{\ x [x,x + he^ x G Q{a 1 h) an elementary 
rectangle of the grid. It is clear by the previous definitions that a vector field / 
is discrete conservative on Q(a, h) if and only if for the counterclockwise directed 
boundary 7 = dE of any elementary rectangle, the discrete curl of /: 

(curhJXz) := ±T^f,j' )h (14) 

= ^ {/ V> * 2 ) + / V + h, x 2 ) + fix 1 + h, x 2 ) 

+ fix 1 + h, x 2 + h)- fix 1 + h, x 2 + h) 

- fix\x 2 + h)- f ix\x 2 + h) -f(x\x 2 )} 

is zero. Observe that a discrete curl is a trapezoidal sum over the edges of an 
elementary rectangle, divided by h 2 , the area of the rectangle. 

Starting with a scalar field g e C 3 (R 2 ), we introduce the following modifi- 
cation algorithm to obtain a discrete conservative vector field / on C*(a, h). By 
translation, we may assume that a — 0. First we set f^x 1 ^ 2 ) = [Djg)ix 1 , x 2 ) 
whenever x 1 = or x 2 = 0. [Dj denotes partial differentiation with respect to 
x j .) 

Let us consider now elementary rectangles of the grid in the first quadrant. 
We proceed inductively with layers of rectangles whose lower left (SW) vertex 
is (x 1 , x 2 ), x 1 A x 2 = rh, r = 0, 1, 2, . . . . Because of symmetry, it is enough to 
describe the algorithm when x 2 < x 1 . In the rth layer we proceed as x 1 = jh, 
j = 0, 1, 2, . . . . The next rectangle inherits the values of / defined on the vertices 
of previous rectangles, except for the upper right (NE) vertex, which is called 
the "new" vertex. For this, compute the modified discrete curl 

(curl»/)(aO := ^ {fix\x 2 ) + fix 1 + h, x 2 ) + fix 1 + h, x 2 ) 

+ iD 2 g) ix 1 + h, x 2 + h) - (D ig ) ix 1 + h, x 2 + h) 
- fix\x 2 + h)- fix\x 2 +h)- fix\x 2 )} , 

(which is not zero in general) and for j = 1 , 2 set 

fix 1 +h,x 2 + h) = iD 3 g)(x 1 + h,x 2 + h) + i-iy^h (cax\ 9 h f)(x) (15) 

at the "new" vertex. It is clear that the so defined / is discrete conservative in 
the first quadrant. Observe that the two modification terms (— l) J_1 /i(curl|/)(x) 
have the same absolute value, and in (curl/j/)(x) they both have minus sign. 

In other quadrants the situation is analogous to the case of the first quadrant, 
but the "new" vertex is the upper left in the second, the lower left in the third, 
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and the lower right in the fourth quadrant. The signs of the modification terms 
±h (cm\fj)(x) in (fTij]) have to be changed accordingly as well. 

It remains to see how large the difference between / and V<? is. For the sake 
of simplicity, consider only points in the first quadrant, points in other quadrant 
being analogous. We claim that these errors accumulate only diagonally. 

In fact, if one considers two neighbor rectangles with a common edge, then 
we can see from (TH)) and (fTS|) that the discrete curl of the "new" rectangle 
(which is right or up from the "old") does not inherit the modification terms 
±h (curl^/) (x) of the "old" rectangle. The reason is that the modification terms 
of the "old" rectangle cancel in the "new" curl; out of the two pairs of edges 
joining at the NE vertex of the "old" rectangle, one pair of parallel edges has 
opposite directions, so the sign of the modification term changes, while the other 
pair has identical directions, so the sign of the modification term remains the 
same. On the other hand, in the case of two rectangles with a single common 
vertex (so which are in diagonal position), the curl of the "new" (NE) rectangle 
does inherit the modification terms of the "old" (SW) rectangle, because both 
pairs of parallel edges have opposite directions. 

Thus (fl5|) implies that for any n > 1 and R > 0, 

sup \f(x)-(D jg )(x)\ (16) 

|a;i|A|x 2 Hn/i;|^|<_R 

R 

< nh sup |(curlfoVsO(x)| < — = sup \(cur\ h Vg)(x)\. 

\x\<R v2 \x\<R 

Thus the error estimation reduces to an estimate between the "true" curlVg 
:= -Di2<? — D%\g — and the discrete curl^ Vg. Or, more precisely, between the 
"true" path integral /^(Vg, "f' )ds = and the discrete path integral T 7 (Vg, 7q)/i 
of the conservative vector field V ' g over the boundary 7 of an elementary rect- 
angle. 

Now, as it is well-known, if cf) € C 2 (R), the error between the integral and 
the trapezoidal area of <fi on [x, x + h] is 

* +h xt ^ Mx) + (t>{x + h) h a .,„ 

cb{u)du~h^ p '- = -—<t>"{x + sh), (17) 

where < s < 1. 

When E := [x, x + hex] x [x, x + he^ is an elementary rectangle and 7 = dE, 
it follows that for any g g C 3 (R 2 ) and for any i£l 2 , 



(curl h Vff)(aO = -^T^Vg,^)h 



1 

1? 
1 

J? 



T 7 (Vg,Y )h- / (V 5 ,7^)ds 



7 



= ^ { {Dlg^x 1 + ath, x 2 ) + (D^g)(x 1 + h,x 2 + s 2 h) 
- {DlgKx 1 + s 3 h, x 2 + h)- (Dlg)(x\x 2 + s 4 h)} , 

where < Sj < 1. This implies that 

|(curlhVs)(a)|<^e B) (18) 
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where e E := sup eE { \ (Dfg) (x)-(Dlg)(y)\,\ {Dig) {x)-{Dlg){y)\}, which goes 
to zero as h — > + . 

Combining (fTB|) and ([T5| . we obtain the following 

Lemma 3. Lef a € R 2 , h > 0, i? > 0, and g £ C 3 (R 2 ). Define e(h) = 
e g {h,a, R) by 

e{h) := su V {\(D]g)(x)-{D]g){y)\ : \x-a\, \y-a\ < R+h; \x-y\ < hV2;j = 1,2} 

(19) 

(which goes to zero as h — > + ). Let f denote the discrete conservative vec- 
tor field on the grid Q(a,h) = a + hi?, obtained from V<? by the modification 
algorithm described above. Then 

sup \f(x)-(Vg)(x)\ < §he(h). 

\x-a\<R 

This lemma expresses the fact that the error of the above modification al- 
gorithm even when divided by h can be made uniformly arbitrary small on any 
bounded planar set by choosing a small enough h. 

5 Planar Ito's formula as an almost sure limit 
of the discrete formula 

Let us apply now the planar (d — 2) case of the discrete Ito's formula ([§]) to a 
random, time-dependent scalar field g : ft x M. + x M 2 — > M, g(u>,t,x), which is 
measurable in w for all (t,x), and is C 1 ' 3 in (t,x) for almost all u>. 

More exactly, fixing a £ R 2 and taking m = 0,1,2,... and h = 2~ m , 
construct first a sequence of discrete conservative vector fields f m {^, t,x) on the 
grids G{a, 2~ m ) from (V<?)(w, t, x), for each uj and t fixed, by the modification 
algorithm discussed in the previous section. (In this paper Vg = (Dig, D^g) 
denotes gradient of g{uj,t,x) with respect to x and D t g its derivative with 
respect to t.) Fixing a bounded time interval [0, if], by a slight generalization 
of Lemma [31 for any uj fixed we get that 

sup sup \f m (uj,t,x)-(Vg)(uj,t,x)\ < f 2~ m e° K (2' m ), (20) 

te[0,K] \x-a\<R O 

where 

e° K (h) sup{\(D 3 jg )(uj,t,x) - (D?g)(uj,t,y)\} 

and the supremum is taken for all | x — a |, \y— a \ < R+h, \x — y\ < h-\/2, t £ [0,K] 
and j = 1, 2. Then e% (2" m ) ->0asro->oo. 

Moreover, it is clear from the modification algorithm (|T5|) that at any point 
of the grid, f m differs from V ' g by a finite linear combination of D^g values, so 
f m (uj,t,x) is continuously diffcrcntiablc as a function of t, like V<?. Moreover, 
by taking derivative of (|15p with respect to t, we obtain 

[Dtf^u^x 1 + h,x 2 + h) = [DjDtg)^^^ 1 + h,x 2 + h) 

+ (-ly-'h (curl^(A/ m ))(x). 
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By the same argument that lead to Lemma [31 for any fixed uj it follows that 

R 

sup sup \{D t f m )(uj,t,x)~(D t \7g)(uj,t,x)\<-he 1 K (h), (21) 

te[0,K]\x-a\<R 

where 

ejf(ft) :=sup{|(^ 3 Ag)(w,i,x) - (Dfl> t 0)(w,f,l/)|} 

and the supremum is taken for all \x— a\, \y— a\ < R+h,\x—y\ <hy2,te [0,K] 
and j = 1, 2. Then e x K (h) -> as h -> 0+. 

Second, start with a planar Brownian motion (W(t))t^M + constructed as 
in Section [21 but shifted so that W(0) = a. Then take the planar Skorohod 
embedded random walks (B m {t)) tm+ , B J m (r2- 2m ) = W j (s 3 m (r)) (j = 1,2) in 
©. That is, let S r := B m (r2" 2m ) and 

X r = X m (r) := 2 m {B m (r2" 2m ) - B m ((r - l)2- 2m )} . 

Then (X m (r))^ 1 is a two-dimensional, independent, (±1,±1) symmetric coin 
tossing sequence. Define stochastic sums by 

n 

(/ m (w,u,W0 • W)T := £{/ m (w,*r-i,Bi,(tr-i),^(*r-i)) 2- l Xi(r) 

r=l 

+/^(w,t r _i,Bi,(t P ),^,(t r _i)) 2- m X^(r)}, (22) 

where t r := r2~ 2m and n := [t2 2m \ . (Of course, B m , X m , and W all depend on 
oj, but this dependence is not shown here and below, to simplify the notation.) 
Now the discrete Ito's formula © can be written as 

{f m {u,t n ,x),i {x)) 2- m (23) 

n 

= E T *=* tr) W"^' *r, X) ~ f m (oJ, tr- 1)X )} , To 2-" 1 
r=l 

+ (/ m ( W , M ,i¥).^)r 

2"^(r) 

/m MmAM) ~ /m (^^i--l,^m(^T-),^mftr-l)) 

2-™Xl{r) 

Our strategy is that we show that each term in this formula, except for the 
stochastic sum, almost surely uniformly converges to the corresponding term of 
the planar Ito's formula, on any bounded time interval. Then it follows that 
the stochastic sum almost surely uniformly converges as well (to the stochastic 
integral), on any bounded time interval. Hence at the same time we obtain a 
proof of the planar Ito formula as an almost sure uniform limit of the discrete 
formula. 

Theorem 1. Suppose g(u,t,x) is measurable in u> for all (t,x), and is C 13 
in (t, x) for almost every lo. For each m — 0, 1, 2, . . . and for all u and t, let 
f m (uj,t,x) denote the discrete conservative modification of (V g)(u>,t, x) on the 




Self-intersection local time based on random walks 



13 



grid Q{a,2 m ). Taking a planar Brownian motion W, for each m define the 
Skorohod embedded random walk B m . Then for arbitrary K > 0, 



sup 

te[o,K] 



(f m (u,u,W)-W) T t 



(Vg)(u,u,W(u)) ■ dW(u) 







almost surely as m — ?• oo, and for any t > we obtain the planar Ito's formula 
as an almost sure limit of the discrete formula 



g(GJ,t,W(t))-g(u,0,W(0)) 



(D t g)(uj,u,W{u))du 



(24) 



(Vg)(u,u,W(u)) ■ dW{u) 



(Ag)(oj,u,W(u))du. 



Proof. We are going to prove (j24|) pathwise. For this, let fio, IP{^o} = 1, denote 
a subset of the sample space f2, on which, asm->oo, B m uniformly converges 
to W on [0, K] and g(to, t, x) is C 1 ' 3 as a function of (t, x). During the proof we 
fix an uj 6 Qq. Then, obviously, W has a continuous path and its range over 
[0,K] lies in a ball Br(o) := {x : \x — a\ < R} with a finite radius R = R(to). 
Also, by Corollary []] we may assume that the range of B m over [0, K] lies in 
the same ball for any m > mo(w). 

Consider the term on the left side of (|2"3")l . We want to show that it uniformly 
converges to g(cv,t,W(t)) — g{u,t,a) for t € [0, K]. Define the path j n = 
[a,(B^(^),a 2 )] + [(B^(t n ),a 2 ),B m (t n )}, where n := [t2 2m \ and t n - n2~ 2m . 
Then we have 



sup 

te[o,K] 



T^Hfm(u,t n ,x), 1 ' Q (x))2- 



W(t) 



< sup \T ln (f m (w,t n ,x),^' (x))2 
te[o,K] 



((Vff)(«,i,x),->i(rs)) ds 
T Jn {{Vg){cj,t,x), 7 ' (x))2-" 



+ sup 

te[o,if] 



sup 

te[o,KT] 



r 7 „((V ff )( W ,t ja ;), 7 ^x))2- m - / ((V ff )(w,t, i),^^)) da 

((V ff )(w, *,!), y (»)) ds . (25) 

Bm(t„) 

2/)/ 



The first term on the right side of (1231) can be bounded above by K2 
2~ m x (i?/6)2- m e? f (2- ro ) = O(e K (2- rn )). Here the first factor bounds the 
number of terms in the trapezoidal sum, the second factor is a multiplier in 
each term, and the third factor bounds the difference of the terms in the two 
sums by (|2U|) . 

The second term on the right side of (|23|) can be bounded by K2 2m x 
(M 3 /12)2- 3m = 0(2- rn ). Here the first factor bounds the number of terms 
in the trapezoidal sum and the second factor bounds the difference of a trape- 
zoidal term and the corresponding integral by (|17l) , where M3 is an upper bound 
of the magnitude of third x-partial derivatives of g for (t, x) £ [0, K] x Bn(a). 

The third term on the right side of (|25|) can be bounded by Mi x(0(m 1! 2~ : %-) 
+2~ m ) = 0(m32~Tr). Here M 2 is an upper bound of \Vg\ for (t,x) € [0,K] x 
Bji(a), while the first term in the parentheses bounds |W(i) — B m (t)\ by ©, 
and the second term bounds \B m (t) — B m (t n )\. 
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In sum, the the term on the left side of (gSJ) is bounded by O (e° K (2~ m )) + 
0(m%2~^). 

Let us turn now to the first term on the right side of (|23p . Define the paths 
j r = [a, (B m (t r ),a 2 )} + [(B m (t r ),a 2 ),B m (t r )} (1 < r < n). Then we have 

n 

E T ^a tr) {{fm(u, tr, x) - f m ( U , t r ^,x)} , ^(x)) 2~ m 
n 

= E T 7r *r-i + 4, x),D t f m (u, tr-i + s 2 ,x)} 2- 2m , y (a;)) 2-" 1 

n 

= E T ^ ({P*I>iff(w,*r-i + 8i,x),DtIhg{uJ,tr-i + s^x)},^)) 2 - m 2- 2m 

r=l 

+ 0(e 1 K (2- m )) 

71 

= E A {Ty„ (Vg(cj,t r ,x), 7 ' (x)) 2~ m ) 2- 2m + 0{e l K (2~ m )) + 0(2- 2m ) 

n „ 

= E D * / V -9( w ' *) ds 2_2m + °(4(2" m )) + 0(2-" 1 ) 

r=l Jfr 

n 

= E A tr, Bm(*r)) - U, a)} 2~ 2m + 0(^.(2^)) + 0(2^) 

r=l 

= / (D t g)(uj,u,W(u))du-g(uj,t,a)+g(uj,0,a)+O(e 1 K (2- m ))+O(mi2-^). 
Jo 

Above we made use of the fact that all considered functions are uniformly contin- 
uous over the bounded set [0, K] x B^(a). The first equality used the mean value 
theorem in the time variable, component-wise for f m , with < s*, s 2 < 2~ 2m . 
The second equality applied inequality (|21j) . combined with the largest possi- 
ble number of terms in the sums and the corresponding multipliers. The third 
equality estimated the error, when one replaces the values of the components of 
\7g at time t r _i + s J r by their values at time t r . The fourth equality replaced the 
trapezoidal sum by the corresponding integral, using (|17[) . In the fifth equality 
we evaluated the integral over the path j r . The last equality replaces the sum 
of the function D t g in the time variable by an integral, and at the same time 
replaces B m by W, using ([5]). 

Finally, let us consider the last term in (|2U|) . By (|2U|) , for the first term in 
the braces we have 

fm ^ tr-l, B m (t r ), gj^tr-l)) ~ f m (u, -Bm(*r-l))) 

2-m X m (r) 

(Dig) (^,^ r _ 1 ,^ft r ),^ft r _ 1 )) - (gig) {m 7 t r _ l7 B m (t r ^)) 

2~™X m (r) 

+ O(e K (2- m )) 

= (Dug) (w, t r -i,Bi(tr-i) + sl,B 2 m (t r ^)) + 0(e° K (2- m )), 

where |s*| < 2~ m . For the other part in the braces of the last term of (j2"3")l 
involving one can obtain a similar result by the help of Z?22- In sum, with 
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|4| < 2~ m , we have 

/m ( W >*r-1,-B m (ir),-B^(£ r _i)) - f m ((J, t r _ 1( B m (t r _l)) 



E 



2"™X m (r) 

/to (^,t r -i,B m (t r )) - / m (w,t r _i,B^(t r ),B^,(t r _i)) 



-2m 



2- ro ^(r) j 

= ]T {(£>115) (W, tr-l, B^(tr-l) + 4,^m(*r-l)) 
r=l 

+ (£ 225 ) (c,^!,^^),^^!) + 4)} 2 - 2 - + 0(e° K (2" m )) 
t 

(Ag)(o;,w,W^(u))du + 0(e? { (2- m )) + (9(mi2-T). (26) 



Here, similarly as above, the sum was replaced by an integral and B m by W. 

Thus we have seen that all terms, except for the stochastic sum, converge 
to their counterparts in (124)) . almost surely uniformly on [0, K]. (An extra term 
—g(ui,t,a) has appeared too on both sides, that cancel each other.) Therefore 
the stochastic sum must converge to the stochastic integral in the same sense 
as well. This ends the proof of the theorem. 

□ 

The reader may have noticed in the statement of Theorem Q] that the usual 
condition in Ito's formulae that the random function g(u, t, x) be adapted to the 
filtration of Brownian motion W, was not needed: the assumed smoothness of g 
together with the pathwise, integration by parts stochastic integration technique 
made this assumption unnecessary. 



6 Discrete self-intersection local time 

The definition of discrete self-intersection local time follows the lines of the 
definition of discrete local time ©, (fTTj) and (|12p . Take first a planar simple, 
symmetric random walk (5 l ra )^ =0 i = 0, S n = Y^i=i -^i ( n — 1)> where (S^) 
and (S%) are independent one-dimensional simple, symmetric random walks 
with unit steps in unit time. Define self-intersection local time of the random 
walk by 

ai(n,x) := : < i < j < n, Sj - <Sj = x} 

TL— 1 TL— 1 

= EEwi. ( 27 ) 

z—0 j=i 

where n £ Z + and i£Z 2 . We also need partial self-intersection local times 

a%(n,x) := : < i < j < n, Sj - Si = x, Sj+i - Sj = fj,}, (28) 

where fj. € { — 1, l} 2 - 

Clearly, by the strong Markovian property of random walks, each inner sum 
£i(n — i,x) :— Y^jZi ^-{Sj~Si=x} m the last term of (|2"T)l is a local time of a 



16 



T. Szabados 



random walk started from the point Si, taken at time n — i at the point x; 

n-l 



Q!i(n, x) = £j(n — i, x). (29) 

i=0 

Denote the largest number of visits to a point of the random walk in the first n 
steps by 

t(n) := sup £(n, x). 
zez 2 

Similarly, denote the largest number of visits to a point of the random walk 
starting from point Si, in the first n — i steps, by t* (n — i). Then for any uj G fi 
one clearly has 

*o(n)>*I(n-l) >■••>£-!(!)• (30) 

In a classical paper [B] , Erdos and Taylor showed the following inequality for 
the maximum number of visits of a random walk in the first n steps: 



n-fr-oo*" log 2 (n) ~~ 7T 



lim sup-; — < — a.s. (31) 



The next lemma is an easy consequence of this result. 
Lemma 4. 

sup xgZ 2 ai(n,a;) 1 
hmsup 2 5 — — a - s - 

n->-oo n log 71 TT 

Proof. By (EH), (1201) and (EH), 

sup„ (: y2 ai(n, x) ,. ?t£n(n) 1 

lim sup Pa:gz 2 < lim sup — < - a.s. 

n-Hx fl log n n->oo n log n I" 



□ 



I do not know if this lemma is sharp or what sharp lower limit rate could be 
given for sup 2 , eZ 2 cti(n, x). In this regard it can be mentioned that Dembo et 
al. [5] relatively recently proved the conjecture of Erdos and Taylor that in fact 

l im = I a . s . 

7WOO log (n) 7T 

Now we apply the previous results to shrunken random walks. Let h > 0, 
x £ hi? = Q(0,h) and t G h 2r L + . Consider a simple, symmetric random walk 
(SVi)^Lo' ^0 = 0, on the grid Q(0,h), with time steps h 2 . (That is, the time 
between step n and step n + 1 of the walk is h 2 .) Define the corresponding 
self-intersection local time as 

a h {t,x) := h 2 : < i < j < t/h 2 , S 3 - S, = x} 

n— 1 n — 1 n— 1 

= h2 E E Hst-s^} =h 2 Y^^- h x/h), (32) 

i=0 j—i i—0 
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where n = t/h 2 and £i (n — i,x/h) is defined in the same way as above. A partial 
self-intersection local time in the direction /i G {1,— l} 2 is 

n— 1 n— 1 

a" h (t,x) :=h 2 ^Y^l {s ._ s . =x . s . +1 _ s . =llh} . (33) 

i—Q j—i 

Let us extend ah(t,x) for any t G R+ and x G R 2 as a continuous function. 
First, with i G h 2 Z, + fixed, we apply linear interpolation in x. Let i be a 
point in a lower triangle A with vertices (a 1 , a 2 ), (a 1 + h,a 2 ), and (a 1 , a 2 + h) 
for some a G hi?. Let A = a/j(f, (a 1 , a 2 )), -B = a/i(t, (a 1 + /i, a 2 )), and C = 
ah(t, (a 1 , a 2 + h)). Then define 

11 2 2 

a h (t,x):=A+^^-(B-A) + ^^-(C-A). (34) 

Analogous is the case with an upper triangle. 

Second, define at{t,x) :— ah(h[t/h\,x) for t G R+ and x G K 2 . Similar is 
the extension of partial self-intersection local times as continuous functions. It 
will be of use later that then 

r h 2 

a h (t,x)dx = —(A + B + C). (35) 
J a 6 

Lemma |4] clearly implies that 

,. sup xeR 2 a h (t,x) 1 
hmsup =-x < — a.s. 

h^0+ t bg 2 (t/h 2 ) ~ 7T 

Briefly, this means that 

sup a h (t, x)=0 (\og 2 {h)) (h -> 0+) a.s. (36) 

xeK 2 

Take now a planar Brownian motion {W(t))t>o, W(0) = 0. Then take planar 
Skorohod embedded random walks (B m (t)) t > ~Bl l (r2- 2m ) = W j {s j m {r)) (j = 
1, 2) for m G Z+. Clearly, B m (t) is a shrunken random walk with h = 2~ m . For 
sake of simplicity, let us denote the corresponding self-intersection local time by 
a m (t,x), and partial self-intersection local times by al^(t,x). It follows that 

svtp a m (t,x) = 0(m 2 ) (m — ¥ oo) a.s., (37) 

uniformly on any bounded time interval t € [0, K]. 

It will be also useful in the sequel that for a.e. u, a m (t,x) — if |x| > 
R = R{uj), for any m > and t G [0,if], supposing R is large enough. This 
follows from the fact that for a.e. uj G f2, the continuous function W(v) — W(u) 
is bounded on the compact triangle Vk = {(u, v) : < u < v < K} and 
B m (v) — B m (u) almost surely uniformly converges to it on Vk as m — > oo. 

7 A discrete Tanaka— Rosen— Yor formula 

The aim of this section is to give a discrete version of the planar Tanaka- Rosen- 
Yor formula. Beyond its intrinsic interest, a special case of this formula will serve 
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as a basic tool for a rather natural definition of planar Brownian self-intersection 
local time. Like the discrete Ito's formulae in Section [3J this formula will be an 
algebraic-analytic one with no intrinsic randomness involved, though, naturally, 
it will be applied in a probabilistic context afterward. 

Lemma 5. Let <fi be a C 3 scalar field in the plane. Fix an h > and let x £ 
hid = Q(0,h). Consider a sequence X r = (X^,X 2 ) (r > 1), where X'l = ±1. 
Take partial sums Sq = 0, S n = h(X± + • • • + X n ) (n > 1), supposing that the 
steps of this "walk" are performed in time units h 2 . Let us take the discrete 
paths 7 r = [0, 0)] + [(S^, 0), S r ], (1 < r < n). Then with any y € hi? fixed, 
one obtains the following discrete Tanaka-Rosen-Yor formula; 

n 

£{T 7b CVftx-Sj-vMoWh} h 2 

n 

= E< 7 V (v<Kx-Sr-vW (x))h} h 2 

n r— 1 

+ EE (^-i - Si - y) hXl + {D 2 4>) ((S^S 2 ^) - Sj - y) hX 2 } h 2 

r=l j=0 

1 n r— 1 

+ - E E( A ^ ( fi r-i ~ Sj -y)h i + 0(h) + 0(e(h)), (38) 

r=l j=0 

where e(h) — > as h — > 0. One has the following equality for the last term as 
well: 

n r — 1 

L h <j>(t n , y ) : = e E( A< w o^-i - s i - y) hi 

r=l j=Q 

= E <Xh{t n ,x)(A<t>)(x-y)h 2 + 0(hlog 2 h), (39) 

where t n — nh 2 and ah(t n ,x) is the self-intersection local time \2T§ of the sums 
S n . For any K > fixed, the error terms in 138\) and 139\) are uniform while 
t n G [0,K]. 

Proof. Define the following time dependent scalar field g v : h 2 l + x hi? — » K, 

t/h 2 

g y(t,x) :=J2^-Sj-y)h 2 , 

3=0 

where y € hi? is a parameter. Take a finite R > such that the disc Br(0) 
cover all the points (S r ) r — 0™, n = t/h 2 and the point y as well. Then all 
points Sj — Si — y are contained by the disc B3r(0). Thus by Lemma [3J one 
can construct a discrete conservative vector field ip in the plane such that 

sup Mx) - V#c)| < f he(h), (40) 

\x\<3R z 

where e(h) = e^(h, R) — > as h —> 0. 
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Further, define 

t/h 2 

fy(t,x) :=Y,^-S 3 ~y)h 2 . 

j=o 

Then by (|40|) it follows that for any y and t fixed, 

sup \fy(t,x)-v g y(t,x)\ 

\x\<3R 

= sup 

|x|<3ii 

and f v (t,x) is a discrete conservative vector field in the plane. 

Now apply the discrete Ito's formula ([9]) to f v (t,x). Let us denote t r = rh 2 
(r = 0, 1, . . . , n). Then we get the following terms. 

The term on the left side of © becomes 

Tts (P(t n ,x) n ' (x))h 

n 

= E^ - Si - y),7o(^)) M ^ + O(eW)- 

i=o 

The error term is obtained since there are t/h 2 terms in the trapezoidal sum, 
there is a multiplier h in each term, and we can apply (|41[) to each term. 
The first term on the right side of (O becomes 

n 

E T *=So (OT-Z) - /^-l,*)} ,7o(*)) ^ 
r=l 

n 

= E ^ - Sr - y),7oW) h) h 2 + 0{e{h)). 

r=l 

The error term is obtained since there are t/h 2 terms in both sums, respectively; 
there are multipliers h and h 2 in each term, respectively; and we can apply (|40j) 
for each term. 

The second term on the right side of © becomes 

n 

EKf) 1 (tr-l,S r -l)hXi + iff (i^!, ^ r 2 } 
r=l 

n r— 1 

= E E U 1 ^) (S-i - Sj - 2/) «^ + (Atf) fi?_i) - S, - y) hX 2 } h 2 

r=l j=0 

+ 0(e(h)). 

Again, the error term is obtained since there are at most t/h 2 terms in both 
sums, respectively; there are multipliers h and h 2 in each term, respectively; 
and we can apply (00]) for each term. 



t/h 2 

E^ a 

3=0 



t/h 2 

Sj -y) ^ 2 -E v< ^ ~ s i ~ y) h " 



3=0 



<^he{h), (41) 
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Finally, the last term on the right side of © of ©becomes 

' (fy) 1 (t r _ 1; - (P) 1 (tr-l,Sr-l) 



E 



(P) 2 (i r _!, S r ) - (.p) 2 (i r _ X , S 2 _ x )) 



hX 2 

7i r — 1 



EE 

r=l j"=0 



(13x0) ((S*, S?^) - Sj -y)- (gig) (S r _i - Sj - y) 



hXf j 

Here the error term is obtained because there are at most t/h 2 terms in both 
sums, respectively; there are multipliers h 2 in both, respectively; each term is 
divided by h; and we can apply P0"|) for each term. 

We need to write the last term in two different ways. The first way mimics 
the method applied to the last term in the proof of Theorem [TJ There exist 
si, s 2 E [—h, h) such that the last term equals 

n r— 1 

£E {{D xl <t>) {{SU + si, S 2 ^) - Sj - y)) 
r=l j=0 

+ (D 2 24>) ((S?,S?_i + si) - S : -y)}h 4 + 0(e(h)) 

11 r— 1 

= J2J2( A ^ ~ Sj -y) h 4 + 0(h) + 0(e(h)). (43) 

r=l j=0 

Here the error term 0(h) is obtained when one replaces the translations s J r by 
in the second partial derivatives of <f>, which are uniformly continuous over the 
bounded ball £?3A.(0). 

The second way of writing the last term (|4"2")l uses discrete self-intersection 
local times (1331) . arranging the terms in (|42l) according to x = SV-l — Sj: 

E <(tn,x)u^x-y)h 2 +0(e(h)) 

xehZ 2 /^6{-l,l} 2 

= Y, a h (t n ,x){A<p)(x-y)h 2 + 0(hlog 2 h) + 0(e{h)), (44) 

where 

{Dx(j>){x l + h^,x 2 ) - (D 1 ^)(x 1 ,x 2 ) 



u M (x) := 



h/j, 1 

(D 2 c^)(x 1 + hfx\x 2 + hp?) - (D 2 (t>)(x 1 + hn\x 2 ) 



h[i 2 

= (DuMx 1 + s 1 , x 2 ) + (D 22 (t>)(x 1 + h^,x 2 + s 2 ) = (A4>){x) + 0(h), (45) 

s 1 ,s 2 g [—h, h}. In (|4"4"|) the error term 0(h log 2 h) is obtained when one replaces 
the translations s 3 r by in the second partial derivatives of 4> in , harnessing 
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the upper bound (j3"6"|) for ah{t ni x) and the fact that there are at most t„/h 2 
non-zero terms (multiplied by h 2 ) in the summation for x. 

Then, collecting the terms of the discrete Ito's formula, by (|43| and (144)) we 
obtain the two versions claimed in the lemma. □ 

It has to be emphasized that so far in this section all obtained formulae 
have been algebraic-analytic ones, independent of any randomness. Now take 
a planar Brownian motion W{t) and replace the sums above by a Skorohod 
imbedded sequence: S n = B m (n2~ 2m )), with h = 2~ m . Then taking limits of 
the discrete Tanaka-Rosen-Yor formula as m — ¥ oo, one obtains a continuous 
version of the formula [20 ;i Theoreme 1] . 

Theorem 2. Suppose that <fi is a C 3 scalar field, W(t) is a Brownian motion 
in the plane, W(0) = 0, and y £ M. 2 . With m = 0, 1, . . . and h = 2~ m , apply the 
discrete Tanaka-Rosen- Yor formula i38\) to the imbedded random walks S n — 
B m (n2~ 2m )). Then, as m — > oo, each term of the discrete formula almost surely 
tends to the corresponding term of the following continuous formula, uniformly 
on any bounded interval t £ [0, K]: 

t 

<f»(W(t)-W(u)-y) du 

ft rv 

= t<p(-y) + (V0) (W(v) - W(u) -y) du - dW{v) 



o 



JO 



t 

(A<f>) (W(v) - W{u) - y) du dv. (46) 



The last term can be written as an almost sure limit of sums involving self- 
intersection local times of imbedded random walks: 

L<j>(t,y):= [ [ (A4>) (W(v) - W(u) - y) dudv 



o 



= lim V a m (t,x) {A^){x-y)2' 2m . (47) 

711 — ^ no ' * 



m— »oo 

x£2 



Proof. To prove the almost sure convergence of the terms in (|46[l , it is essentially 
enough to apply Theorem [1] The only new element here is that in each term 
there is a Riemann sum of a continuous function, that converges to the corre- 
sponding Ricmann integral for almost every path as m — > oo. The value of the 
parameter y £M. 2 should also be approximated by a closest point y m £ 2~ m I? , 
for which \y — y m \ < 2~ m . This does not cause any problem, since all functions 
in (|46p are continuous in y. 

Thus the limit of the term on the left side of (1381) is 



t 

{<() (W(t) - W(u) -y)-<j> (-W[u) - y)} du. 



n 



The limit of the first term on the right side is 
t 

{</> (W(u) - W{u) -y)-<f> (-W(u) - y)} du. 
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The extra J Q — (— W(u) — y) du term appears on both sides, so can be can- 
celed. Observe that / Q * — T4^(u) — y) = t4>{—y). 

The limits of the second and the last terms of the right side of ([55]) are 
clearly the corresponding ones in (|46l) . The equality in (|47p clearly follows from 
the equality ^ . □ 



Since <j) G C 3 (IR 2 ) and is a.s. continuous, it follows that, almost surely, 
the term on the left side and the first and the last terms on the right side of (|46p 
are continuous functions of (t, y) G R+ x R 2 . This implies the same conclusion 
for the second, stochastic integral term as well. 



8 A definition of planar self-intersection local 
time 

A possible definition of ordinary local time in one spatial dimension uses a 
special case of Tanaka's formula applied with the function <j){x) — x V 0, which 
is a fundamental solution of the one-dimensional Laplacian d 2 /dx 2 ; see this kind 
of definition for example in [TS1 p. 117] and [3J Section 7.2]. The definition of 
self-intersection local time presented below is a suitable planar modification of 
it. This means that our definition uses a special case of planar Tanaka-Rosen- 
Yor formula with the function <p(x) = log|x|, which is a fundamental solution 
of the planar Laplacian A = Du + D22, ignoring a constant multiplier. 

Let W(t) be a planar Brownian motion, W(0) = 0. For each x € M 2 and 
5 > 0, we define an everywhere continuously differentiable approximation of 
log \x\ by 



4> S (x) := 



Mj^ + log^ iov\x\<6, 
log [as j for |a;| > 6. 



Then 



and 



(A<f> s )(x) 



■§! for \x\ < 6, 
for \x\ > 6. 



Note that (A<^> )(x) is not defined for \x\ = S, but we set it to be there. 

Since <j) S is not C 3 , Theorem [5] is not directly applicable to it. However, by 
a standard procedure, taking a convolution with a sequence of C°° functions 
q n with compact support shrinking to {0}, and then taking a limit as n — > 00, 
solves this problem. For sake of explicitness, let q(z) = cexp (—(1 — M 2 ) -1 ) f° r 
|z| < 1 and otherwise, where the constant c is chosen so that J R2 l{ z ) dz = 1. 
Put q n {z) = n 2 q(nz) and <f) S n = (f> s * q n (n > 1). 

Then <f> 5 n G C°°(R 2 ); <j) S n -> cp 5 , V<p s n -> V0 5 both uniformly in R 2 ; while 
A4>n — > A(j) S pointwise except for [as j = 5. Thus one can apply Theorem [5] to 
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and take a limit of the terms as n — > oo. The resulting formula is 



^ - W(u) - y) du 



o 



= t(j)\y)+ ( [ {V(f> s )(W(v)-W(u)-y) du- dW(v) 
Jo Jo 

+ 1 A ({(u, ») e y t : - W(u) - 2/| < <5}) , (49) 

where A denotes planar Lebesgue measure, Vt = {(u, v) : < « < v < t}, 
and y € M 2 . It also follows that each term here is an almost surely continuous 
function of (t, y): this is clear for each term except for the second, stochastic 
integral term on the right side, but then it follows for this term too. 
It is important that, by (147[) . the last term can be written as 

\l4\t,y) = ^ lirn^ ]T a m (t,x)2- 2m , (50) 

xe2- m Z 2 nB 6 (y) 

where Bs(y) is the closed disc centered at y with radius S. 

Rosen |12] suggested the following definition of self-intersection local time of 
planar Brownian motion W. Define the following occupation measure for plane 
Borel sets A and time t > 0: 

fj, t (A) :=\({(u,v) : < u < v < t,W(v) - W(u) e A}), 

where A is planar Lebesgue measure. Rosen [TJ] then proved that the self- 
intersection local time a(t, x) := a.s. exists when x ^ 0. 

An alternative approach is to consider the symmetric derivative of /it w.r.t. 

A: 

/. \ ,. Vt{B&{y)) 

a(t,y) = hm — — — 

v :i " 5^o+ X(B 5 (y)) 

= lim — A (Uu,v) : < u < v < t, \WM - W(u) - y\ < 5}) (51) 

where B$(y) denotes the disc centered at y with radius 5. Among other things, 
the next theorem establishes the a.s. existence of this finite symmetric derivative 
for any t > and y =t 0. It is well-known (see e.g. Rudin [14]) that when there 
exists a finite symmetric derivative of fit w.r.t. A except for the point 0, then 
Ht is absolutely continuous w.r.t. A on M 2 \ {0} and the symmetric derivative 
equals the Radon-Nikodym derivative; the support of the singular part of /i t 
can only be the point 0. 

Theorem 3. The terms of \4ty almost surely converge to the corresponding 
terms of the following Tanaka-Rosen-Yor formula as S — > + ; when y ^ 0: 

t 

log|W(i) - W(u) - y\du 

tl0gM+ / / \W(v)-W(u)-y\ 2 dU - dW{v) + * ^ (52) 



o 
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cf. \2(A (2.j)]. Moreover, all terms, including a(t,y), are a.s. continuous in 
(t, y) when y 7^ 0. 

It also follows that the self-intersection local time a(t,y) of planar Brown- 
ian motion is the almost sure limit of averages of self-intersection local times 
ctm{t,y) of imbedded random walks: 



a(t,y)= lim lim V a m (t,x)2 2m 

xe2- m Z 2 C\B s {y) 

= lim lim — — / a m (t,x)dx {y^O). (53) 

<5_>0+ m^oo 7TO J Bs (y) 

Proof. The well-known properties of planar Brownian motion imply that for any 
t > and y ^ 0, inf{|W(£) - W(u) - y\ : < u < t} > 0, with probability 1. 
Hence it follows the almost sure convergence of the left side of (|4^|) as 6 — > + 
when y ^ 0. Moreover, the integrand converges monotonically as S — > + , so 
the limit and the integral can be interchanged. It also follows that the left side 
of (|52"|) is a continuous function of (t, y) when y ^ 0. 

The convergence and the continuity of the first term on the right side is 
trivial when y ^ 0. . 

By Lemma [7] in the Appendix, the second term on the right side of (|49l) is 



o J o 



(V0 4 ) (W(v) - W{u) ~y)du- dW(v) 

1 f* f v f W[v)-W{u)-z 



Jo Jo Jb s ( v) \W(v) - W{u) - zP 

if ( r* r v w{v) - w{u) - z 



dzdu ■ dW(v) 
- du ■ dW(v) \ dz. (54) 



™* 2 JB s (y) Uo Jo \W(v)-W(u)-z\ 

The interchange of integrations is allowed by an extension of Fubini theorem. 
(Remember that by Lemma [5j the stochastic integral is an almost sure limit 
of discrete sums.) By Lemma |5] in the Appendix the stochastic integral has a 
continuous version for y ^ 0. Thus by the mean value theorem of integrals, 
([51)1 has an almost sure limit if y ^ as 8 — > + , namely the one stated in 
the theorem. By Lemma [Etc), the limit is a continuous function of (t,y) when 

The above limits imply that the last term on the right side of (j4U)) has an 
almost sure limit na(t, y) if y ^ 0. (l53l) follows from this by ([50]) and Lemma 
[5] in the Appendix. It also follows from the above arguments that a(t,y) is 
continuous in (i, y) when y ^ 0. 

□ 

Formula ([53"[l is the definition of planar self-intersection local time which has 
been the main objective of the present paper. It is an open question if the limits 
in (|53[) can be interchanged; then one would get the more impressive almost sure 
limit a(t, y) = lim^oo a m (t, y). 

The most interesting question is 'What happens to the planar self-intersection 
local time a(t,y) when y —> 0?' It was discovered by Varadhan in 1968 that a 
goes to oo then, and, in fact, it has a logarithmic singularity at 0. So one can 
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introduce renormalized self-intersection local time 7 by the formula 
ot(t,y)- ^log-^r when y ^ 0, 



\y\ 

lim^o a(t, y) - ± log ^ when y = 0. 



It was shown by Le Gall [TU] that the limit above exists both almost surely and in 
L 2 . Below we will give an alternative proof of it, together with a Tanaka-Rosen- 
Yor formula for 7, based on Theorem [3J Yor [2D] (2.k)] proved convergence in 
probability by a similar approach. 

Prior to that, let us give the expectation of 7, cf. |10) . 

Corollary 2. 

*y(t,v) = ( ^^-^H-^)-^ Wh6n V * ' (55) 
{ £(k>g(2t)-C-l) when y = 0, 

where Ei denotes the exponential integral function and C is Euler's constant. 
Thus this expectation is finite and continuous for every (t,y) £ K+ x K 2 . 

Proof. In formula (|52")l the expectation of the left side is given by Lemma EJa) , 
while the expectation of the stochastic integral term is by Lemma Ufa), see 
the Appendix. □ 

Since j(t,y) has finite expectation for any (t,y), quite often in the litera- 
ture the renormalized self-intersection local time is defined by subtracting its 
expected value. Since that would complicate some formulae below, here we do 
not follow that practice. 



Theorem 4. Combine the first and the last terms on the right side of t52\) into 
a 7 term. When y — > 0, the resulting terms in A5ty) converge almost surely and 
in L 2 to 



t i-t 



, log\W(t)-W(u)\du= I I V ™1 V) ™1 U) du- dW(v)+n 7 (t, 0). (56) 
Jo Jo \W(v) - W{u)\ 2 

Moreover, "/(t,y) is finite-valued and continuous in (t,y) even when y = 0. 

Proof. Let us denote the stochastic integral in (l52l by Y(t, y), cf. ((67)) in the 
Appendix. Define Y n (t) :— Y(t, (n^ 1 ,0)), n > 1. By Lemmata), for each n, 
Y n {t) is a continuous L 2 -martingale. By for any k > 1, 

E\Y 2k (t)-Y 2 „-i(t)\ 3 <c Q {K)k 5 2- 3k (0<t<K), (57) 

where cq(K) is a constant depending only on K. 

Thus a basic martingale inequality implies for any k > 1 that 

p( sup \Y 2k {t)-Y 2k -i{t)\>2- k /A <2 3k / 2 E\Y 2k {K)-Y 2k -i(K)\ 3 

\0<t<K J 

< c Q {K)k b 2- 3k ' 2 < Cl (K)2- k , 

where c\{K) is a constant depending on K. 
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An application of the Borel-Cantelli lemma then yields that Y 2 k (t) almost 
surely uniformly converges for t £ [0, K] to a continuous L 2 -martingale as k — > 
oo. (The convergence is also in L 2 by (|57|l .) Using isometry, the expression 

11 W» - - (2- fc ,0) , 

du 



|W^)-iy( U )-(2-fc,0)| 2 

converges in L 2 as well when k oo. Since for any fixed v, W(v) — W(u) 
almost surely does not equal to when < u < v, here the integrand can be 
dominated for any large enough k > fco( w )- So the integral and the limit can be 
interchanged and we get that the limit of Y 2 k is the continuous L 2 -martingale 
Y(t, 0), which is the first term on the right side of (|56p. 

Moreover, Y(t, y) converges to Y(t, 0) as well when y 0, a.s. uniformly for 
t € [0, if]. For, by an argument based on (|6"9")l . similar to the above one, given 
any e > 0, for any y, 2~ k ~ 1 < \y\ < 2~ k with k large enough, 

\Y(t,y)-Y(t,0)\ < \Y(t,y)-Y 2k - 1 (t)\ + \Y 2k - 1 (t)-Y(t,0)\<e 

a.s. uniformly for t € [0, if]. 

The convergence as y — > of the term X(t,y) on the left side of (j52|) can be 
treated analogously by Lemma Hfb) of the Appendix. □ 

The following occupation time formulae, cf. |10| . follow from the previous 
results. 

Corollary 3. Suppose that f : M 2 — > R is a bounded, Borel measurable function. 
Then 

t 

f(W(v)-W(u))dudv= I f(x)a(t,x)dx 



o Jo 



/(.D^^^-ilogH^ dr. ir„si 



Alternatively, 



o Jo 



{/(W(«) - - E/(W(«) - W(u))} dudv 

f{x) {j(t, x) - Ej(t, x)} dx. (59) 



Proof. It is enough to show (|58|) for indicator functions of discs; by standard 
methods, linear combinations of indicators extend to a general /. So let us take 
a closed disc B r (a) and show that 

t r-v i> 

/ l Br(a) (W(v)-W{u))dudv= / a(t,y)dy 

JO JB r (a) 



[ j{t,y)dy-- [ log\y\dy. (60) 

JB r (a) K JB r (a) 
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From the results above we know that the stochastic integral term on the right 
side of (1521 has a version Y(t, z) which is continuous for any (t, y) £ [0, K] x R 2 . 
([55]) gives 



t rv \ f 

(V/) (W{v) - W(u) -y)du- dW{v) = — / Y(t, z) dz. 
o Jo nd Jb s ( v ) 

Substitute this into (gH]): 

* i r 

<f> s (W(t) - W(u) -y)du- tcj, 6 (y) - — / y(t, z) dz 



i 

^2 



l Ba( ,,)(W»-W»)dud<;. 



'0 JO 

Then integrate this equality over a closed disc B r (a) with respect to y: 

' / (W(t) - W(u) -y)du- t(p\y) - -L [ Y(t, z) dz \ dy 

B r (a) \J0 ™ JB 5 {y) I 



1 

J 2 " 



/ / l^ te) (W(«)-W( U ))dydudu. (61) 

JQ JB r (a) 



Now, it is clear that 

If 1 

r 1 ™. - ?2/ 1 B^(y)( z ) d J/ = (a) (*) + o^B.fajfz), (62) 
<5^0+ 7T0 JB r (a) 1 



where B°(a) = {x : \x — y\ < r}. Take limit in (fBTj) as S — > + . By Thcorem[3l 
the limit of the expression in the braces is na(t, y) when y 7^ 0. Since a(t, y) = 
7(i, y) — ^ log |y|, and 7 is a.s. continuous everywhere while J B ^ log \y\ dy is 
finite even if £ B r (a), this gives us the right hand sides of (j60|) . 

The left hand side of (ITOf is obtained as the limit of the right hand side of 
(I5T|) , using (f6"2"j). since it has zero probability that W(v) — W(u) £ dB r (a). This 
completes the proof of (|55|) . 

Now take expectation in ([55 



E/(W(w) - W(u)) dudv= I f(x) <j E 7 (t, x)-- log |x| \ dx. 



'0 ./0 

Subtract this from ([55| . and the result is (15T)1) . 



7T 



□ 



Finally, we add some remarks. We can combine the occupation time formula 
55|) and formula (f^)]) . with y = 0, when is a C 3 scalar field: 



(W(t) - W{u)) du- / (V</>) - J¥(u)) du • dW(v) 

Jo Jo 



(A(f>)(x)a(t,x) dx+t(l)(0) 



1 
2 

i/ (A^)(a;)(7(t,x)--log|i|^ ds + t#0). (63) 

2 ./ D 2 7T 
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Green's theorem implies that (f)(0) — j- L 2 (A(j))(x) log \x\ dx holds when — 
x ■ (S7(fj)(x) log \x\ tends to as \x\ — > oo. In this case (|6"3"j) simplifies to the 
following Tanaka-Rosen-Yor formula: 

<f) (W(t) - W(u)) du= / (V0) (W(u) - W(u)) du • dW» 



^0 JO 



(A(j))(x)"f(t,x) dx. 



Further, comparing the first equality of ([53"f with (|4"T1) results the almost 
sure limit 

lim V a m (t,x)(A<j))(x)2- 2m = [ a(t,x) (A<f>)[x) dx. (64) 

Now, almost surely, the support of the continuous a, and then by ©, the 
support of each a m for to > toq(w), can be covered by a finite disc -Br(O) with 
a large enough radius i? = R(u). Thus Lemma [6] implies that the sum here can 
be replaced by an integral. Also, the Poisson equation A(f> = f can be solved 
in the plane for any continuous /, f(x) = 0(\x\~ 2 ~ 5 ), 5 > 0, by the formula 
4>(y) = ^- J R2 (A(f>)(x) log \x — y\ dx. Thus, for any such /, (jMf can be written 
as 

lim / a m (t, x) f(x) dx = I a(t,x)f(x)dx a.s. (65) 

m-yoo J R2 J R2 

This weak convergence formula supplements the basic definition (I53[) of self- 
intersection local time. 



9 Appendix: Some technical lemmas 

This lemma says that the sum in (|50[) can be approximated by an integral. 
Lemma 6. Almost surely, for any t € R+, j/£l 2 and 5 > fixed, 

a m (t,x) T 2m = f a m {t,x)dx + 5 0(m 2 2~ m ). 

Proof. First, let us estimate the error between the sum of a discrete function 
and the integral of an interpolated function over a rectangular domain A n — 
\xQ,Xn) x [yo,y n ) with vertices on a grid hi 2 , in general. So let / : hi 2 — > R 
be a discrete function and f(x,y), (x,y) e M 2 , be obtained from / by linear 
interpolation on triangles of the grid, as it was described for ah in (|34D . Put 

n n 

S n -=h 2 ^ f{ x hVi) = tf^y^fixi-i^j-x), 

{xi, yj )eA„ i=l 3 = 1 

and 

T n ■= / f(x,y)dxdy 

JA n 

h 2 " n 

= ~Q~ ^2 + 2 /( X *> + 2 f( x i-liVj) + /(^i: %)} > 

i=l j=l 
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cf. (|55|) . Then in the error T n — S n , all contributions of inner vertices cancel 
and only the contribution of vertices at the boundary of A n remain: 

h 2 " -1 

T n - S n = — ^ {f(Xi,Vn) + f(x„,yi) ~ f(xo,Vi) ~ f{Xi,Uo)} 
i=l 

h 2 

+ y {2/(x , Vn) + 2/(x n , y ) + f{x n ,y n ) - 5/(x , y )} ■ 

Returning to the statement of the lemma, by our definition in Section H2 
a m (t, x) is obtained by linear interpolation on triangles of the grid with mesh 
h = 2~ m . If one considers the difference of the sum and the integral over a disc 
B$(y), the cancelation of inner vertices still holds, and only the contributions 
of vertices adjacent to the circumference remain. The number of these latter 
vertices is of the order of 6 0(2 m ). Thus by (|37p we have that 



^ a m (t,x)2 2m - / a m (t,x)dx 



< 6 0(2 m ) 0(m 2 ) 2- 2m = S 0{m 2 2- m ). 
This completes the proof of the lemma. □ 

We need the following representation of 
Lemma 7. 

(V^)(x) = 4a / r~\ dz 
nS 2 Jb 5 (o) \x-z\ 2 

= — \z f ]og\z-x\ n (z) ds(z) (iel 2 ), (66) 



Ci(0) 



where Bs(0) is the closed disc centered at the origin with radius 5 > 0, Cg(0) is 
its counterclockwise directed boundary, hq is the outward unit normal along the 
boundary, and ds denotes integration with respect to arc length. 

Proof. The second equality follows from a standard theorem of vector analysis; 
note that the discontinuity of the integrand in the second term when \x\ < S is 
not essential. To prove that the first term equals the third, because of rotational 
symmetry, it is enough to consider points x = (— a, 0), a > 0. Then the third 
term of (|6"6")) becomes 

~^sj 1 °g(( a + <5cos0 ) 2 + ( <5sin0 ) 2 ) (cos#,sin#) d9= f-^* (f ) . 



where 



u for < u < 1, 
for u > 1. 



i r ( 

: = — / Iog(w 2 + l + 2ucos6l)cos6ld(9= <^ " 
These prove the equality with V(f> given by (HHJ). □ 



30 T. Szabados 

The next lemma establishes some important properties of the stochastic 
integral appearing in Theorem [3] 

Lemma 8. Fix an arbitrary K > 0. Consider the stochastic integral 
VU ^ f* f V W(v)-W(u)-y 

Y(t > V) = Jo Jo \W(v)-W(u)-y\> du ' dW{v) - (67) 
Then the following properties hold. 

(a) Y(t,y) is a continuous L 2 -martingale with expectation as a function of 
t e [0, K] for any fixed y el 2 , 

0>) 

E\Y(t,y)-Y(t,y>)\ 3 <C\y-y'\ 2+ P (68) 

with a finite C = C(K, a) and with an arbitrary f3 G (0, 1) for any t € [0, K] 
and \y\, \y'\ > a, where a > is arbitrary, fixed. More exactly, 

E\Y(t,y) - Y(t, y')\ 3 < c(K) (l + log + 2KOT 2 ) log 4 ^—j-^y - yf, (69) 

\y-y\ 

where c(K) is a finite constant depending only on K and log + x := OVlogx. 

(c) Y(t,y) has a version which is a.s. a continuous function ofy^O. In fact, 
it has a version which is a.s. a continuous function of (t,y) when y =/= 0. 

Proof, (a) Since for any fixed v, W(u) := W(v) — W(v — u) is planar Brownian 
motion as well that starts from 0, we have 



E\Y(t,y)\ 2 = / E 







W(v) - W(u) - y , 
■ du 



o \W(v)-W(u)-y\ 



2 



dv 



t 1 - - ••< - 2 



< / E / -= du dv = tE[ - — — r du 

~ Jo \Jo \W(u)-y\ J \J \W(u)-y\ 

Then by symmetry and by the independence of increments of W, we get 
that 

t r t ^ 



E\Y(t,y)\* <tE[ dm du 2 



o Jo 



■\W( Ul )-y\\W(u 2 )-y\ 



2t / dui dzi : r / dii2 dz2 



[0,t]xB 2 ^ui\zi- y\ J [ui t]xR 2 ' 2tt(u 2 - ui)\z 2 - y\ 

Writing z 2 = y + r(cos 0, sin#), z\ = y + p(cos a, sin a) and v — u 2 ~ u\, for 
the inner integral here we obtain 



f*-" 1 dv r 27r de r°° dr 

Jo v Jo 27r Jo 



J2 , 2 



+ p* -2-rp cos (fi 



dv r* de f°° dr (r- P ) 2 







1 y/vJo ^ Jo V2^ - V V U 
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Thus 

/.2tt 



E\Y(t,y)\ 2 <2t / / / dpe-£r2V27r(i- Ul ) 

< 8?r^ 2 < oo (0 < t < K, y E R 2 ). 

(b) First, by the Burkholder-Davis-Gundy inequality, for any m > 0, there 
exists a finite c m such that 

/ \ I m 



E|Y(t,y)-y(t,2/)| 



* - W(it) - y W{v) - W{u) - y' 



o \W(v) - W{u) - y\ 2 \W(v) - W(u) - y'\'- 



du 



dv 



Using the elementary vector equality \a\a\ 2 — b\b\ 2 | = \a — 6|(|o||6|) 1 , 
and the fact that for any fixed v, W(u) :— W(v) — W(v — u) is planar 
Brownian motion as well that starts from 0, it follows that 

E\Y(t,y)-Y(t,y')\ m 



<c m \y-y'\ m ^ \ ((( ^ du] dv\ 

\Jo \Jo \W(u)-y\\W(u)-y'\ J J 

= c m t^\y-y'rE( f —— 1— -du) , (70) 

\Jo \W(u)-y\\W{u)-y'\ J 

where for any fixed v, W{u) := W(v) — W(v — u) is planar Brownian motion 
as well that starts from 0. 

Thus to show (|68p . it is enough to give a suitable upper estimate for the 
last expectation in (|70[) when the norms of y and y' are bounded below by a 
and t £ [0, K]. Now, by symmetry and by the independence of increments 
of Brownian motion, with m = 3 we obtain that 

•* 1 N 3 

h ' / \W(u)-y\\W(u)-y'\ dU 



pt pt ft 3 

6E / dui / du 2 / du 3 J | 

JO J Ui J U2 j — i 



LL \W(u 3 )-y\\W(u 3 )- y >\ 



duidzi— / dit2dz2- 



[0,t]xl 



27ra i J[u u t]x®.* 27t(u 2 -wi) 

l«3 I 2 3 
e -2(„ 3 - U2 ) J* X 

du 3 dz 3 — -II- -. (71) 

u 2 ,t] X R2 2tt(m 3 - u 2 ) f^j \Zj - - y'\ 

Without loss of generality, from now on we may assume that < \y — y'\ < 
1/2. B r (x) will denote the closed disc centered at x with radius r. Here 
and later we use the following covering: 

N 

R 2 = (S 1 (y)Ui? 1 (y')) C U (J (C n (y) U C n (y')) , 
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where C„(y) := Ba^-w^v) n S 2»-i|„-„'| (f) n B c 2n - 1]y _ yl] {y'), 0<n<N, 
N = flog(|y-y|-7log2l, and C^(y) := B 2 -^_ yl \{y). For C n (y') the 
definitions are similar. 

To show the method, let us estimate the innermost integral I± in (1711) using 
the above covering of R 2 . First, 



112 



I{ := / du 3 / dz 3 



g 2( U3 - U2 ) 1 



(Bi^uBife')) 27r(w 3 - u 2 ) \z 3 - y\\z 3 - y'\ 



t — U 2 



2 



<l dv dz 3 '—— = t-u 2 . 

Jo Jri 2ttv 

Second, write z 3 = y + r (cos 9, sin 9), z 2 = y + p(cosa,sina), and for n = 
0,1, ... ,N obtain that 

i , 3 —i I 2 

' e 2 <"3-«2) 1 



Ji,n(2/) := / du 3 / dz 3 



"2 



c„ 27r ('"3 - u 2 ) \z 3 - y\\z 3 - y'\ 



2 - d9_ f'-^dv dr e~ r +P - 2 Z- iB - a) 



o 2tt J v J 2n -i\y-y>\ ' 2™ 1 |y — y'l 



2^ 



d0 f f - U2 dv /•(2 n |»-*'|-/'costy, v- , - 

: e a* / ds- 



o 2vr 7 y^; y(2~- 1 |i,- V '|- / cicoBe)/v^ 2 ™ l \y-y'\ 

2k d0 ,t-« a d« e _ £ _W9 



2tt Jo w 

(The value of a clearly does not matter, so it was replaced by 0.) Here one 
can use the simple estimate 

1 _j> 1 x 

-e « dv < - + log, — . 
v e o 

Then 

2(t-u 2 )\ d0 



h,n{y)<J Q (~+l°s 



+ p 2 sin 2 9 J 2tt 



1 2(t- M2 ) f 2 \ . 2 d9 



< 2 + log + 1 

(n = 0, 1, . . . , N), since p = \z 2 — y\. The estimate for Ii- 1 (y) is the same, 
and so is for any I\, n (y') replacing y by y' . 

In sum, the estimate for the innermost integral is 

JV 



<t-u 2 + 12\og—^— f 2 + log + 2{t U2) ), (72) 

\y-y \ V \%2-y\\z2-v\J 
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since N + 2 < 61og(l/|y - y'\) when < \y - y'\ < 1/2. 

The estimation of the second and third integrals in (|7ip can go in a similar 
fashion. Omitting the details, the result is 



3 



E U \W(u)-y\\W(u)-y'\ dU . 

1 / It 
< C (X)log 4 , l + l og+ — 

\y-y'\ V \v\\y 

< c{K) log 4 (1 + log + 2Ka- 2 ) , (73) 

\y -y\ 

for any y,y' such \y\, \y'\ > a, < \y — y'\ < 1/2, and t e [0, JC], where 
c(if) is a finite constant depending on K. By (|70)) this verifies with 
arbitrary j3 < 1, and so proves the lemma. 

(c) By (b), Y(t,y) as a function of y satisfies the condition of a special case of 
the Kolmogorov-Chentsov theorem, so a.s. it has a continuous version as a 
function of y when y =/= 0. 

One can similarly show that Y(t,y) has a version which is a continuous 
function of (i, y) when y ^ 0, e.g. considering fourth moment instead of the 
third. 

□ 

This last lemma investigates the properties of the integral on the left side of 
formula ([52 



Lemma 9. Consider the integral 

X(t, y) := [ log \W(t) - W[u) -y\du (t>0,y€ M 2 ). 



Then it has the following properties, 
(a) 

EX( t ,y) = tlo g \y\-^Ei(-^)-he-^ {y + ), (74) 



2t J 2 

limEX(t,y) = EX(t,0) = -(log(2t) - C - 1), (75) 

y^O 2 

where Ei denotes the exponential integral function and C is Euler 's constant. 

(b) 

E\X(t } y)-X(t,y')\ 3 <C\y-yf + P (76) 

with a finite C = C(K, a) and with an arbitrary j3 € (0, 1) for any t € [0, K] 
and \y\, \y'\ > a, where a > is arbitrary, fixed. More exactly, 

E\X(t, y) - X(t, y')| 3 < c(K) (l + log + 2Ka' 2 ) log 4 J^q/\\y ~ y'f, (77) 

where c(K) is a finite constant depending only on K . 
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Proof, (a) For any fixed t, W(u) := W(t) — W(t — u) is planar Brownian motion 
as well that starts from 0. Thus 



EX(t,y) = / log\W(u)-y\du 



log \x - y\ 
2iru 



e 2 « dxdu 



r f r r 3 f d6 1 
= / du dr-e~^ / — - log (r 2 + p 2 - 2rpcos(9 - a)) , 
Jo Jo u Jo 2lT 2 

where a; = r(cos 6*, sin 6) and y = p(cosa,sina). It is clear that the last 
integral does not depend on a, so we can replace a by 0. Since 



d6 1 
2^2 



log (r 2 + p 2 — 2rpcos(9)) = log(r V p), 



it follows that 



EX(t,y) 



du 



dr—e 2 « log(r V p), 



that gives exactly the results (|T4")) and (|75|) . 



(b) 



E|X(t,y)-X(t,y')| 3 =E 



log 1 - ' — — du 



o \W(u)-y'\ 



< E 



lot 



\W(u)-y\ 



\W{u)-y'\ 



du 



Using the inequality | log& — loga| < (a A — a| < (a -1 + — <A 

for a, & > 0, and then symmetry and the independent increments of W , we 
obtain that 



'\|3 



l\X(t,y)-X(t,y')\ 



< \y-y'\ 3 ® 



(iW^-yl^ + lWiu)-^]- 1 ) du 



<6\y-y'\ 3 / duidzx- 

i[0,t]xK 2 /7TU 1 J[u lt t]> 



du 2 dz 2 



e 2(u 3 -u 2 ) 



g 2( U2 - U1 ) 

27r(w 2 - «l) 



1 



Since this last formula is very similar to formula (|7ip . the remaining part 
of the proof is essentially the same, so omitted. The result differs from the 
one of Lemma [H](b) only by a constant multiplier depending on K. 

□ 
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